Nonlinear Wigner Solid Transport Over Superfluid Helium Under AC Conditions 
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Nonlinear transport properties of the two-dimensional Wigner solid of surface electrons on super- 
fluid helium are studied for alternating current conditions. For time-averaged quantities like Fourier 
coefficients, the field- velocity characteristics are shown to be qualitatively different as compared to 
that found in the DC theory. For a spatially uniform current we found a general solution for the 
field-velocity relationship which appears to be strongly dependent on the current frequency. If the 
current frequency is much lower than the ripplon damping parameter, the Bragg-Cherenkov reso- 
nances which appear at high enough drift velocities acquire a distinctive saw-tooth shape with long 
right-side tails independent of small damping. For current frequencies which are close or higher 
than the ripplon damping coefficient, the interference of ripplons excited at different time intervals 
results in a new oscillatory (in drift velocity) regime of Bragg-Cherenkov scattering. 
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INTRODUCTION 

Electrons trapped on the surface of liquid helium form 
a clean two-dimensional electron system of which the av- 
erage Coulomb interaction energy can greatly exceed the 
average kinetic energy (for a review, see [l|, 0]). For 
typical electron densities realized experimentally (n s < 
10 9 cm -2 ), the Fermi-energy of surface electrons (SEs) 
is much smaller than temperature. Therefore at a low 
enough temperature, depending on the surface electron 
density (T c cx -^/nj), this electron system undergoes a 
transition to the Wigner solid state. This was first ob- 
served by Grimes and Adams Q from the onset of res- 
onances induced by electron interaction with capillary- 
waves (ripplons) whose wave- vector q is close to electron 
reciprocal lattice vectors g. Since the electron lattice 
determines a specific set of frequencies ui n = uJ Sn [here 
n = 1, 2, 3, ... , uj q = a/ a/pq 3 / 2 is the ripplon spectrum, 
a is the surface tension, and p is the mass density of liquid 
helium], experimental evidence for a triangular electron 
lattice on a liquid-He surface was given. Electron inter- 
action with such ripplons appeared to be very strong, 
leading to a huge reconstruction of the WS phonon spec- 
trum in the low frequency range 

The resonances of Grimes and Adams occur when the 
frequency of the input signal lu is close to u> n . Transport 
properties of the WS of SEs are usually studied under 
much lower frequencies (lu <C loi). Nevertheless, even 
under low frequency conditions the resonant interaction 
with ripplons of frequencies which are close to w„appear 
to be possible as a nonlinear conductivity effect [5(. The 
physics of this phenomenon can be explained as follows. 
The pressure at the interface induced by the WS mov- 
ing with a constant drift velocity v can be represented 
as a series of terms proportional to exp [zg- (r — vt)} , 
where r is the in-plane coordinate vector. As a func- 
tion of t, it can be considered as a series of harmonic 



perturbations with frequencies gv, and one can expect 
a resonance response of the system when gv is close to 
u>„. It should be noted that the corresponding velocity 
vi = uJi/gi is rather low (typically about 1 — lOm/s ) 
while the thermal velocity of electrons in the liquid state 
is usually much higher (about 3 • 10 3 m/s). Therefore the 
single-electron Cherenkov emission of ripplons is a quite 
usual phenomenon contributing into resistivity of SEs on 
liquid helium. The important point is that for a moving 
WS with gv —>L) n , the response of the system can be con- 
sidered as a coherent Bragg-Cherenkov (B-C) scattering. 
This effect was first described by the usual perturbation 
treatment Q which leads to symmetrical peaks of the 
electron collision rate with non-Lorentzian tails. 

Besides direct B-C scattering effects limiting the WS 
velocity Q , there are other interesting nonlinear conduc- 
tivity (a) phenomena observed for SE transport over su- 
perfluid helium. In the presence of the magnetic field 
oriented normally to the surface, function of 

the input voltage V has a remarkable N-type anomaly 0] • 
The decreasing part of cr~^(V) was attributed to the B-C 
scattering. When studying WS confined in the channel 
geometry, periodic conductance oscillations with vary- 
ing the drift velocity were observed [8|. These oscilla- 
tions were attributed to anisotropic spatial order with 
lines of electrons along the channel edges. Complicated 
nonlinear conductivity of the WS was observed for cur- 
rent frequencies which are close to typical frequencies of 
plasmon- ripplon coupled modes 0]. Interesting WS ve- 
locity jumps caused by the decoupling of electrons from 
the surface deformation were recently observed for the 
WS in a channel [1(3] ■ 

Unfortunately, for q —>■ g\ the electron-ripplon cou- 
pling is strong which leads to a huge increase of the elec- 
tron effective mass at low frequencies due to surface dim- 
ples [4( ■ Under such conditions the perturbation treat- 
ment is doubtful and coupled WS phonon-ripplon modes 
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are usually treated in a self-consistent way [llj, ll2j . In 



this treatment the most of the interaction Hamiltonian of 
the WS with ripplons of q — > gi is included in the descrip- 
tion of the coupled modes. Therefore, a simple classical 
model of coherent B-C scattering of capillary waves for 
the WS moving with a constant velocity [l3| seems to be 
more appropriate than a perturbation treatment. This 
model was introduced in order to analyze a complicated 
nonlinear magnetoconductivity observed previously Q. 
An extension of this model applied to liquid 3 He allows to 
explain the nonlinear field-velocity characteristics of the 
WS under strong ripplon damping conditions 14, 1 15 ) . 

Application of the models of coherent B-C scattering to 
the nonlinear WS transport on liquid 4 He is difficult for 
several reasons. First, experimental geometries usually 
imply that the driving electric field and electron current 
density are not spatially uniform. Secondly, measure- 
ments are done under AC conditions when the electron 
velocity and driving electric field are periodic functions of 
time with the period 2?t /to which is much longer than the 
typical ripplon oscillation period 2ir/u>i. Moreover, the 
damping of ripplons 7 g in superfluid 4 He is anomalously 
small [16] . which means that B-C peaks of the classical 
model are extremely narrow and some other effects not 
included in the model can significantly affect its main 
results. 

In this work we study the surface-displacement profile 
and field-velocity relationship for spatially uniform alter- 
nating motion of the WS over superfluid 3 He and 4 He 
under small ripplon damping conditions: j gi -C u>x- For 
arbitrary frequency of the current to, the exact expression 
for the field- velocity relationship can be found. This solu- 
tion appears to be strongly dependent on the ratio oo/^g. 
Therefore, we separate two frequency regions: lo <C 7 3 
and lo > j g . In both regions, the nonlinear field- velocity 
characteristics obtained here for time averaged quantities 
differ significantly from those obtained in the DC model. 
In the high frequency region 7 9 < to <C lo 9 , we expect the 
appearance of a new B-C scattering regime of the WS 
transport caused by interference of ripplons excited at 
different time intervals. This frequency region is usually 
realized in experiments on nonlinear WS transport over 
superfluid 4 He, and, therefore, we expect that our new 
results will help to understand the nonlinear electronic 
response observed. 



DIMPLE PROFILE EVOLUTION INDUCED BY 
THE DRIFT VELOCITY 



We investigate shape variations induced by the WS ve- 
locity and ripplon damping which are very important for 
understanding the nonlinear WS transport. 

We assume spatially uniform motion of the WS, which 
means that all electron lattice sites have the same dis- 
placement vector s(t) in external fields. In this case the 
electron pressure at the interface induced by the WS 
moving with an arbitrary velocity can be presented in 
the following form: 



(r, t) = n s £ V g exp [ig • (r - s(t))] , (1) 



where V q = V q exp (— q 2 (u 2 ) /4) , the electron-ripplon 
coupling V q depends on the holding electric field E± 
directed normally to the surface and on the wavenum- 
ber q 0, and (u 2 ) - is the mean-square displacement of 
electrons from lattice sites due to fast coupled phonon- 
ripplon modes whose frequencies are limited by Wf > wi. 
Actually, u>{ is the frequency of electron oscillations in 
the potential of a steady dimple. A simple self-consistent 
treatment [T3| gives 
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where (uq) ~ 1.248S/ (2m e Ct \fitn s ) is the mean-square 
displacement due to zero point vibrations, 
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m e is the free electron mass, and c 2 is the transverse 
sound velocity of the electron solid. 

Considering ripplons as an ensemble of surface oscilla- 
tors with a certain damping parameter 7 9 surface dis- 
placements £(r) induced by pressure perturbations of 
Eq. fl} can be found in a quite general form: 



n s gV g 



sin [cjg (t — t 1 )] x 



P -uj g 

x exp [-ig • s(i') + 7 9 (i' - t)] dt' 



(4) 



where Cj n 



Here wc extend the treatment 



given in Ref. [2| into the range of a finite but small ripplon 
damping. For AC conditions s(t) — Sosin(wt). Introduc- 
ing the new variable r = t' — t and assuming lot <C 1 
we represent s(t') as s(t) + v(<)t. In this approximation 
surface displacements can be found as 



The analysis of the classical B-C scattering given in 
Ref. [HI was restricted to a simplified one-dimensional 
DC model. The damping effects were considered in a 
phenomenological way. Here we consider more realistic 
2D model of alternating motion of the WS with a par- 
ticular ripplon damping defined for both 3 He and 4 He. 
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where v = sqlo cos(wt) is the WS velocity. This equation 
represents dimple sublattice moving in-phase with the 
WS. In the limit gv — > 0, Eq. surely gives the well 
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known shape of steady dimples which is independent of 
damping. If gv ~ v g , the shape of dimples is affected by 
the WS velocity and damping. Thus, under the condition 
lot <C 1, which according to Eq. (@J requires w«7 9 , the 
dimple shape changes with time continuously in such a 
way that it is always adjusted to a given velocity v(t). In 
other words, for any fixed value of v(t) the dimple shape 
is the same as that defined by the DC theory with the 
corresponding WS velocity. 

For liquid 3 He, the weak ripplon-damping regime can 
be realized only for superfluid phase at T < 0.3 mK. In 
this case, -f g is determined by ballistic bulk quasiparticle 



scattering from an uneven interface 15j: 
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where is the Fermi momentum of quasiparticles in liq- 
uid 3 He, A is the excitation gap, and f(x) = (e x + I) 1 . 
The ripplon damping of superfluid 3 He decreases with 
cooling at an exponential rate. Still, in experiments on 
WS it can be just reasonably small (7 9 /w s ~ 0.1 or 0.01). 

In contrast, the ripplon damping in superfluid He is 
anomalously small. In the ballistic regime it is given 
by 
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where U4h c is the first sound velocity. For n s = 10 9 cm -2 
and T = 0.5 K, a simple estimate gives r y gi /u> gi ~ 10~ 4 . 
Thus, for WS transport on superfluid 4 He the damping 
coefficient of ripplons is extremely small. It is remarkable 
that this damping coefficient has the same dependence on 
the wave-vector q as that found for the ballistic regime 
of liquid 3 He [Eq. ©]. 

In the reference frame moving along with the WS, the 
dimple profile given by Eq. ([5]) can be evaluated as 



p|Ac(v)| 



-2gv7 ff sin(gr)], (8) 



where D s (v) 



(gv) — i2gv r y g . At low tempera- 



tures, especially for superfluid 3 He, the summation over a 
large number of g is necessary to ensure the convergence 
of the result. Consider r = (x, 0) and assume that the 
driving force — eE is directed along to the x-axis. Then, 
according to Eq. 1(5)), in the absence of damping, dimples 
have a symmetrical shape with regard to electron lat- 
tice sites. A finite damping introduces asymmetry in the 
dimple shape due to the terms proportional to sin^x). 

Shape variations of surface dimples induced by a fi- 
nite WS velocity are very sensitive to orientation of the 
vector v with regard to symmetry axes of the WS. We 
shall consider the following two typical directions. The 




-10- 



x / c 



5- 



(b) 



u = u = 1 




x / a 



FIG. 1: Variations of the dimple sublattice profile £(x,0) in- 
duced by the WS velocity for two typical velocity orientations: 
NN-direction (a), and SN-direction (b). Steady dimples are 
shown by the solid line. The dimensionless velocity increases 
from u = 0.6 to higher values by steps equal 0.1. Calcu- 
lations are performed for superfluid 3 He, n s — 10 s cm -2 , 
E± = 189 V/cm, and T = 0.277 mK ( ji/uii =0.1). 



direction of the vector v which is parallel to the line con- 
necting two nearest neighbors of the electron lattice will 
be called NN-direction. The direction of v which is paral- 
lel to the line connecting second neighbors will be called 
SN-direction. If the ripplon damping coefficient 71 is not 
too small (about 0.1 • the asymmetrical variations of 
the dimple shape start already at substantial shifts from 
the B-C resonance as shown in Fig. Q] (a and b) for two 
typical directions of the velocity vector. Here the dim- 
ple shape of a motionless WS is shown by the solid line. 
The first shape-line (dashed) of moving dimples is cal- 
culated for dimensionless velocity u — v/vi = 0.6 (here 
Vi = u gi /gi). For each next curve shown in this figure, 
the parameter u is increased by steps equal 0.1 except 
for the curve calculated at u = 2/y/if. The last curve is 
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FIG. 2: Variations of the dimple sublattice profile £(x,0) in- 
duced by the WS velocity near the first B-C resonance for 
very small ripplon damping 71/wi = 0.01. The drift veloc- 
ity is oriented along the NN-direction. Calculations are per- 
formed for superfmid 3 He, n B = 10 s cm' 2 , E± = 189V/cm, 
and T = 0.2 mK. 



calculated for the critical velocity u c of the first B-C reso- 
nance, which equals 1 for the SN-direction (Fig.[lJ>), and 
2/1/7F for the NN-direction (Fig. QJi). For the symmetri- 
cal dimple shape, the average force acting on electrons is 
obviously zero. The strong shape asymmetry appears in 
order to transfer the kinetic friction acting on the dimple 
sublattice by the environment to the electron crystal. 

For significantly smaller damping coefficient ^1 = 0.01 • 
cl>i, which is expected for superfluid 3 He at T = 0.2 mK, 
substantial shape changes appear only if the WS velocity 
approaches the first B-C resonance or exceeds it. Dimple 
shape variations which occur near the vicinity of the first 
resonance are shown in Fig. [2] The important point is 
that the velocity-induced displacements found for u = 
u c = 2/t/tt are much larger (about 15 times) than surface 
displacements in the initial surface dimple. This means 
that dynamic decoupling of the WS from surface dimples 
which shall be discussed below for a fixed field condition 
is accompanied by creation of huge displacement waves 
moving in the same direction. 

At u < u c , in spite of huge changes of the dimple profile 
the average position of an electron remains the same be- 
ing fixed to the potential minimum formed by the dimple 
potential and the driving electric field. For experimen- 
tal conditions with a given current, we may consider the 
evolution of surface dimple profile even at u > u c . For 
example, at u = 1.3 surface displacements induced by 
the WS velocity are already substantially reduced, still 
their amplitude is larger than the initial dimple depth. 

It is instructive to consider dimple shape variations in- 



FIG. 3: Variations of the dimple sublattice profile £(:r,0) in- 
duced by the WS velocity away from the first B-C resonance 
condition for 71/wi = 0.01. The drift velocity is along the 
SN-direction. Other conditions are the same as in Fig. 



duced by WS velocity which is away from the first B-C 
resonance condition. The corresponding calculations are 
shown in Fig. [3] As expected, at u = 0.8 < u c the WS 
velocity just increases the dimple depth. Significantly 
faster velocities u ~ 10 cause the opposite effect: the sur- 
face displacements are substantially reduced. For u = 9 
we accidentally reached one of higher B-C resonances, 
causing strong asymmetry in the dimple shape. A small 
detuning up to u — 10 makes the shape nearly symmetri- 
cal which means that the force transferred to the WS by 
dimples is close to zero. It is instructive that in the limit 
u — > 00 surface displacements still remains as shown in 
Fig. [3] by the straight dotted line, but instead of individ- 
ual dimples we have a row of interface valleys oriented in 
the direction of motion. 

Shape variations induced by B-C resonances lead to 
strong changes in the associated mass of an electron dim- 
ple Md which is given by 



P g (^-( g v) 2 ) 2 + (2 7 gv) 2 ' 

where v is directed along the x-axis. For a finite thickness 
of the liquid helium d, each term in the sum should be 
multiplied by coth(gd). Obviously, the B-C resonances 
increase the associated mass of a surface dimple. In the 
limiting case gv — * 00, the associated mass Md(v) disap- 
pears because the dimple lattice is rearranged in a raw 
of valleys (terms with gv = do not contribute in Md 
because of the proportionality factor g%). 
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FIELD-VELOCITY CHARACTERISTICS 

The asymmetry of the dimple shape with regard to 
the average electron position in a lattice site causes a 
force F( D ) acting on the electron solid. In equilibrium 
this force is balanced by an external driving field. By 
definition, F( D ) is the sum of forces acting on each elec- 
tron — J2 e <9Vint/<9r e averaged over electron distribution 
within the dimple [here Vjnt (r) is the electron-ripplon in- 
teraction Hamiltonian whose Fourier transform V q was 
used in Eq. (JTj)]. Electron distribution caused by long 
wave-length fluctuations with tu < lj\ occurs together 
with surface dimples and therefore it should be excluded 
from the averaging. Then, the average of the electron 
density operator n_ q = J2 e ex P (*1 r e) can be found as 

(n_ q ) f = N e cxp [-q 2 (uj) /A + ici- B (t)] (J q , g . (10) 

Using this equation and Eq. Q in the general expression 
— (J2 e dV\ nt /dr e ) [ , the force acting on the WS can be 
written as 

n aV 2 f°° 
(t) = -N e V / sin (Cj.t) e-~<° T x 

g P^g Jo 

xsin{g [s(t) - s(t - t)]} dr. (11) 

For any spatially uniform displacement s(t) given, this 
equation defines the in-plane force induced by the dimple 
sublattice. 

If the time scale of the WS displacement vector s(t) 
is much longer than t ~ l/7g, then g [s(t) — s(t — r)] 
can be approximated by gv(i)r. In this limit F^ D ' has 
the same form at that given by the DC treatment with 
the constant velocity replaced by v(t). In equilibrium, 
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as well as the kinetic friction of the electron lattice 
-N e v e v caused by electron scattering of other 



kinds are balanced by the external force N e eE (here we 
assume that magnetic field is zero). The solution of the 
balance equation can be represented as a field-velocity 
characteristic E(v): 
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(12) 



Using this equation, E(v) can be calculated numerically 
for any given damping coefficient and the collision fre- 
quency v e caused by electron scattering with thermal rip- 
plons, vapor atoms, or even walls if WS is formed in a 
channel geometry. 

For the DC case, Eq. (fT2]) is a two-dimensional ex- 
tension of the classical one-dimensional model of B-C 



scattering reported previously 13( with the real ripplon 
damping parameter and with more accurate electron- 
ripplon coupling. The driving field found from the bal- 
ance equation has sharp maxima in the vicinity of B-C 



given and the WS velocity is adjusted to the field, then 
regions with dE/dv < are unstable. This means that 
for driving fields exceeding the major maximum of E(v) 
the balance of forces is not possible and the WS decouples 
from surface dimples. According to Fig. [21 decoupling of 
the WS is accompanied by creation of huge surface waves 
moving with the group velocity u c . 

In experiments on WS transport, usually it is the cur- 
rent which is given, while the driving field is adjusted to 
the current by electron redistribution which screens the 
external potential variations. This is supported by the 
fact that regions with dE/dv < are experimentally ob- 
served [1, Q] . Therefore, field- velocity characteristics of 
electrons moving ultra-fast with u > u c are very impor- 
tant for understanding the nonlinear WS transport on 
superfluid helium. 



For liquid He with j g / 
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Eq. QJJ applied 



to the DC case would give just a set of extremely sharp 
peaks. At the same time, beyond these peaks E(v) is 
close to zero. For liquid 3 He, the parameter 7 S / u g can 
be much larger than it is for liquid 4 He (even larger than 
unity) and B-C peaks of E(v) can substantially overlap. 
It is worth noting that E(v) depends strongly on the 
velocity vector orientation with regard to the 2D electron 
lattice. For example, in the NN-direction of motion there 
is only one major peak with |g| = g± at u — , while 

for the SN-direction there are two equivalent peaks with 
|g| = g\ at u = 1 and u = 2. 

In contrast to models of B-C scattering discussed pre- 
viously, the real experimental situation has one unavoid- 
able complication: DC measurements are practically im- 
possible and by now all data are obtained under AC con- 
ditions. This means that WS velocity and the driving 
electric field are periodic functions of time with the pe- 
riod 2tt/u>. The frequency of experimental signal is usu- 
ally varies from 10 4 s _1 to about 6 ■ 10 5 s _1 which is much 
lower than the typical ripplon frequency uj\. Therefore, 
it is conventionally expected that the DC model of WS 
transport should give qualitatively correct description of 
data obtained in such experiments. We shall see that this 
is not true for two major reasons. 

First, we note that the condition lo -C u>i is not suf- 
ficient for adiabatic adjustment of surface displacements 
to WS velocity variations. For example at the B-C res- 
onance condition [u — u c ) we have a huge wave which 
follows the WS and for a change to much smaller dis- 
placements a shape variation (without changing the am- 
plitude) is not sufficient. At the B-C resonance the 
capillary wave accumulates great energy which should 
be transferred to the environment in order to make the 
transition from the surface displacements calculated for 
u = u c ~ 1.155 (see Fig.) to the surface displacements 
calculated for u — 1.3 or even for u = 1.1. This means 
that an adiabatic AC extension of the DC model requires 
an additional restriction: the frequency of the current 
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should be much lower than the corresponding ripplon 
damping (co <C 71). This is actually the condition which 
allowed us to transform Eq. (TTT]) into Eq. (fT2|). For WS 
transport over superfluid 3 He this condition is satisfied 
even at j g / co g ~ 1CP 2 . Regarding liquid 4 He, the con- 
dition co <C 71 requires to use an AC frequency which is 
much lower than 10 4 s . 

Secondly, even if the above noted condition is fulfilled, 
Eq. (fT2|) cannot be used directly for plotting field ve- 
locity characteristics. Under AC conditions, it is im- 
portant which quantities are actually measured and pre- 
sented in the field-velocity characteristics. If time aver- 
aged quantities are considered, in the nonlinear regime 
the outcome can be qualitatively different for different 
kinds of averaging. For example, even for harmonic ve- 
locity uq cos cot, the mean-square time averaging of the 
driving field yj (E 2 ) and averaging of the absolute value 
(\E\) give qualitatively different results for the function 
Eiuocoscot) defined by Eq. (|12|l . This could be easily 
proven by considering the right-side tail of the B-C res- 
onance in the limiting case 7 9 — > 0. For averaging \E\, 
in this limit a part of the corresponding integrand can 
be rearranged as the 8 -function [see Eq. (fT2]) ]. and the 
final result will not depend on j g . In contrast, the in- 
tegrand of the quantity {E 2 ) is squared, and, therefore, 
the resonance tail of y/ (E 2 ) increases with reducing j g . 

As a measure of the alternating field one can choose 
the main term of the Fourier series representing E(t): 



to 



7r/u 



E w = — I E(t) cos (cot) dt. 

TT J-it/lj 



(13) 



It is also a kind of time averaging, and similar to (\E\) 
it has a finite resonance tail in the limiting case 7 9 — ► 0. 
In the following we shall consider only quantities E u and 
(\E\) which will be used for presenting the field- velocity 
relationship. 

Assume that the conditions of given current is realized 
and u(t) = uo cos cot. Using the adiabatic treatment dis- 
cussed above, we insert v(t)= viu(t) into Eq. (112|) and 
evaluate the time integral of Eq. (I13|) . It is convenient 
to introduce two integer variables m and n to describe 
the reciprocal lattice vectors g m n = mgW + ng( 2 ) (here 
and g( 2 -* are primitive vectors of this lattice). Then, 
after some algebra, the Fourier transform E u can be rep- 
resented as a function of the velocity amplitude 

Euj(u a ) = Vi ^j- — y -j- Q(p m , n Uo, p m ,n) + 



ea 3 / 2 qi r }l 2 

a m.n ym.n 



H viit . (14) 
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Here we use the following dimensionless notations 
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Q(w,fl) 
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(1 + y 2 - w 2 ) + w 2 (3 2 (1 + y 2 ) 

(16) 

The integral Q(w,fl) can be evaluated in an analytical 
form: 



Q(w,fl) = ^Re 
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, (17) 



where G(0) = fly/1 — /3 2 /4. Remarkably, in the limit 
lg ~ * (fl —> 0) for w > 1 there is a finite asymptote 



Q(w,fl) 



26(w - 1) 
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where 9(x) is the unit step-function. 

For (\E\) as a function of the dimensionless velocity 
amplitude uq, an equation similar to Eq. (|14[) is found. 
The only difference which appears for such averaging is 
that instead of Q(w,fl) we should use another function 
I(w, fl) defined by 



I(w,fl) = ^Re 



arctan 



• s /iu 2 -l+/3 2 /2+iG 



yjw 2 -1+ fl 2 /2 + iG 



(19) 



As expected in the limiting case "f g 
has a finite asymptote 



(fl -»■ 0), it also 



I(w,fl) 



6(w-l) 



y/w 2 



1 



(20) 



which means that the right-side tails of B-C resonances 
are independent of small damping. 

Already from the analysis of the integrals Q(w, fl) and 
I(w, fl) given above one can conclude that B-C resonance 
tails of the AC treatment differ (even qualitatively) from 
that found for the DC models. The left-side tail becomes 
even steeper for both E^(uo) and (\E\), while the right- 
side tails extend far beyond the resonance and are inde- 
pendent of ripplon damping in the limit j g — * 0. This 
behavior is illustrated in Fig. 0] where E u (uo) is plot- 
ted for two typical directions of the WS velocity [NN- 
direction (solid line), and SN-direction (dashed line)] as- 
suming 7 g — > 0. The other parameters are taken for the 
liquid 4 He case. Thus, instead of 5-peaks of the classi- 
cal B-C scattering model here we have saw-tooth shaped 
peaks with long right-side tails. As noted above, field- 
velocity characteristics depend strongly on the direction 
of the WS velocity. It is interesting that for the SN- 
direction there are two major B-C peaks and the second 
one (at uq — > 2) becomes even more prominent than the 
first one because in the AC cause at uq > 2 the velocity 
sweeps through the both resonances. Of course, consid- 
ering the limiting case r ) g /co\ — > we should keep in mind 
that co should be much smaller than j g . Therefore, real 
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replaced by Q(w,(3,u') defined as 




FIG. 4: The first Fourier coefficient of the driving field vs 
the drift velocity amplitude for two typical orientations of the 
velocity: NN-direction (solid line) and SN-direction (dashed 
line). Calculations where performed for superfluid 4 He, n s = 
10 9 cm -2 , T = 0.5 K, assuming the ripplon damping is set to 



ripplon damping should be taken into account for con- 
sistent analysis of field-velocity relationships induced by 
B-C scattering. 

For nonlinear WS transport over superfluid He, the 
ripplon damping parameter is extremely small and the 
condition 7 g < u <C u g is realized in most of known ex- 
periments. This case requires a special treatment because 
large surface displacements excited at a first instance of 
u(t) = u c cannot be relaxed back within the period of 
current oscillations. Therefore a new steady regime will 
be developed for each AC frequency, which can be far 
away from the solutions found above for the condition 
w < 7 9 - 

We have to return back to the exact solution of 
Eq. pip . In this equation now we insert s x (t) — 
(yo/u) sin (ut) and then perform the time averaging de- 
fined by Eq. (|13[) . Here we disregard the unimportant 
correction induced by v e . Then quite generally, -E^vo) 
can be found as 



9xn s gVg f°° . T (ujt\ 

t-^I I sin (u„t) e 79 T cos — 

epu g J v 9 ' \ 2 J 



2g x — sm — 

U \ 2 



(21) 



where Ji(z) is the Bessel function of the first kind. For 
E u as a function of the dimensionless velocity amplitude 
uo, an equation similar to Eq. (|14p can be found. In the 
general case, the integral Q(w,(3) of Eq. (TT4"]) should be 



Q(w, (3, J) = 4 / sin (2y) e~ py cos (u'y) x 
Jo 
w 



xJi 



2 — sin (u'y) dy 



UJ 



(22) 



where 



The dimensionless function 



Q(w, /?, u') describes the shape of a single B-C resonance 
of E u (uq) for arbitrary current frequency and ripplon 
damping (here w oc uq and f3 oc 7 S ). It is easy to see 
that in the limiting case u' — > and — > analyzed 
above, Eq. (|22|) provides the correct asymptote shown in 
Eq. 1(15]). For conditions 7 5 < w <C w g , and .g^vo > 1, the 
argument of the Bessel function entering the integrand of 
Eq. (|2"Tj) attains huge numbers because sin(o7r/2) ~ 1. 
This can lead to remarkable field-velocity relationships 
with side-oscillations which we shall discuss in the fol- 
lowing. 



RESULTS AND DISCUSSIONS 

Consider briefly the velocity-field relationship Eu(u) 
for WS transport over superfluid 3 He at T = 0.25 mK. 
At such a temperature the damping coefficient j gi ~ 
1.54 ■ 10 6 s _1 is much lower than the B-C resonance fre- 
quency u\ and is still much higher than the typical cur- 
rent frequency u used in experiments. The later condi- 
tion makes the adiabatic AC extension of the B-C scat- 
tering model introduced here applicable. The results of 
numerical evaluations of Eq. (Tl4|) are shown in Fig. [5] 
The velocity-field characteristic given by the DC clas- 
sical model is shown by the dotted line. It consists of 
series of B-C peaks with nearly symmetric tails. For the 
first Fourier coefficient as a function of the drift veloc- 
ity amplitude, the AC theory gives lower peaks which 
are asymmetric with regard to the maximum positions. 
It is important that time averaging used in evaluation 
of E u (uq) do not smooth out the B-C resonances com- 
pletely. Another important feature of the AC treatment 
discussed here is the appearance of long right-side tails of 
the B-C resonances. The left-side tails (uq < u c ) become 
even more steeper, because under AC conditions elec- 
trons spend only a little time near the B-C resonance. 

The most interesting experimental results on nonlinear 
WS transport were obtained employing liquid 4 He @, 0]. 
In this case even for T ~ 0.5 K the ripplon damping pa- 
rameter given by Eq. is extremely small. First, we 
assume that u is low enough to make Eq. (|14p applica- 
ble. The comparison of results obtained for the two kinds 
of averaging of the electric field (E u and (\E\)) is given 
in Fig. [6l assuming that the velocity is along the NN- 
direction and v e = 0. Here the sharp peaks of the DC 
model (dotted line) are strongly smoothed by the time 
averaging of the AC model (E u is shown by the solid 
line, and (\E\) - by the dashed line). Additionally the 
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FIG. 5: Field velocity relationship for the DC case (dashed 
line), and the first Fourier coefficient of the driving field E u 
vs the drift velocity amplitude for the AC case (solid line). 
Drift velocity is oriented along the NN-direction. Calculations 
where performed for superfluid 3 He, n s = 10 8 cm -2 , T = 
0.25 mK. 



FIG. 7: The main Fourier coefficient of the driving field vs 
the drift velocity amplitude for different values of the electron 
collision frequency v e : (solid line), 2.4- fO 9 s _1 (dashed line), 
and 7.5- 10 9 s" 1 (short dashed line). Drift velocity is oriented 
along the NN-direction. The DC case results are shown by 
the dotted line. Calculations where performed for superfluid 
4 He, n s = 10 9 cm -2 , T = 0.5 K. 
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FIG. 6: Field velocity relationship for two kinds of aver- 
aging of the alternative driving field: E^ (solid line) and 
7r (\E\) /2 (dashed line). Drift velocity is oriented along the 
SN-direction. Calculations where performed for superfluid 
4 He, n s — I0 9 cm -2 , and T = 0.5 K with the ripplon damp- 
ing parameter defined by Eq. (0. 



maximum values of the B-C peaks are greatly reduced 
as compared to the results calculated for the DC case 
(dotted line). 

The influence of a finite electron collision frequency 
u e due to scattering with thermal ripplons and walls is 



analyzed for the SN-direction and shown in Fig. [7J A 
reasonable estimate for the electron collision frequency 
v e ~ v\ — 2.4- 10 9 s _1 (for chosen n s , E± and T) is found 
considering electron scattering with thermal ripplons in 
the usual way [H and taking into account that at low 
temperatures the average kinetic energy of an electron in 
the WS state differs substantially from T. For experi- 
ments with WS in the channel geometry, v e can be even 
higher because of the WS friction at the channel walls. 
In order to illustrate this effect is Fig. [7J we considered 
also a larger value v e ~ u 2 — 7.5 • 10 9 s _1 . According to 
this figure, the electron collision frequency v e affects the 
both tails of the field velocity characteristic. At the left 
side v e and time averaging of Eq. (|13[) act in the opposite 
ways. The left-side tail (uq < u c ) becomes less steep for a 
finite v e . At the right side v e acts in the same way as the 
time averaging, increasing the right-side tail and making 
Euj(uo) more flatter in the region uq > u c . In general, due 
to the both these effects the field-velocity characteristic 
of the WS acquires a distinctive "nose" shape. 

The numerical calculations presented in Figs. [5] and [JJ 
were done assuming lu <C 7 Sl . In experiments on WS 
transport over superfluid 4 He this condition was not re- 
alized. The influence of the condition lu > j gi on the 
field- velocity characteristics can be understood using the 
general expressions for -E w (v ) and Q(w,0, to 1 ) given in 
Eqs. (f2"Tj) and (f2"2"]l . The main features of B-C scattering 
under AC conditions can be revealed from the dimension- 
less function Q(w, (3,u>') which describes the shape of a 
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single B-C resonance [see Eq. (|14p]. Consider the main 
B-C resonance when we can set w — uq, (3 = 2ji/uji 
and uj 1 — uj/gJi, and for simplicity assume (3 = 0.1. In 
the limiting case u>' <C f3, the function Q(uq,/3, 0) coin- 
cides with Q(uo,f3) obtained in Eq. (fTT|) . For example, 
it is practically impossible to distinguish Q(u , f3, 0.001) 
shown in Fig. [8] by the solid line from Q(uq,(3) given by 
Eq. (fT7|) . As a function of the dimensionless velocity, 
Q(uq, (3, 0.001) has a typical saw-tooth shape discussed 
above. 

Remarkable shape transformations of Q(uq, f3, a/) as a 
function of uq occur when uj' approaches and exceeds the 
value of the parameter which is proportional to ripplon 
damping. A sharp (from the left-side) saw-tooth peak 
of a single B-C resonance is developed into distinctive 
smooth oscillations which [according to Eq. (jT4j)] result 
in similar oscillations of E u (uo). The amplitude and the 
period of oscillations are gradually increase with uj' in 
the range considered. These oscillations represent a new 
regime of B-C scattering of the WS which occur under 
the AC condition, when the current frequency becomes 
comparable or larger than the ripplon damping. 

The period of new conductivity oscillations depends on 
the relation between the frequency of the current uj and 
the frequency of ripplons excited uj g . According to Fig.0 
it increases with the ratio cu/tOg. The ripplon damping 
just increases the amplitude of oscillations. The later is 
illustrated in Fig. [9] where Q(uq, f3,uj') is plotted vs u 
for a fixed value of the ratio uj/u>i = 0.05 and different 
values of the damping parameter. 

There are other important points which follow from 
Figs. [S] and [5] For low damping the first maximum 
of the field-velocity characteristic can be substantially 
larger than the B-C peak value found in the limiting case 
u> -C 7i. Secondly, due to the finite frequency uj, even 
for a very small ripplon damping the left-side tail is not 
steep as it was for ui <C 71. Additionally, the maximum 
position is substantially shifted to higher drift velocities. 

The physics of side-oscillations in the field- velocity re- 
lationship can be explained as follows. If uq > u c , then 
even during a period of current oscillations the WS passes 
through the B-C scattering point four times. The phase 
difference between surface waves of the same q = g ex- 
cited at different times increases with the velocity am- 
plitude. Thus, depending on the velocity amplitude uq 
the excited waves can interfere constructively or destruc- 
tively, which is the reason for the side-oscillations. The 
higher frequency of the current, the larger amplitude is 
necessary to produce the same phase-difference. 

Most of experiments on the nonlinear WS transport 
are performed using a Corbino geometry of which the 
alternating current is spatially nonuniform, and only a 
limited area of the WS can satisfy the B-C scattering 
conditions. This area changes with time because of the 
AC conditions. This experimental situation is very diffi- 
cult to analyze. One may conclude that spatial variations 
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FIG. 8: Shape transformations of the function Q{uq, f3,u)') 
which describes a single B-C resonance under AC condi- 
tions. Calculations are performed for current frequencies: 
ui' = 0.001 (solid line), 0.05 (dashed line), 0.1 (dotted line), 
0.2 (dash-dotted line), 0.3 (dash-dot-dotted line). 
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FIG. 9: Shape transformations of the function Q(uo, (3, u>') 
which describes a single B-C resonance under AC condi- 
tions for a fixed frequency [lu' — 0.05) and different ripplon- 
damping coefficients: f3 — 0.1 (solid line), 0.02 (dashed line), 
and 0.01 (dash-dotted line). 

of the current would additionally smooth out the B-C 
resonances. It should be noted that the field-velocity 
characteristics of the WS with a "nose" shape without 
the B-C peaks where observed in the experiment on WS 
transport in a channel geometry 0] . It is interesting that 
conductance oscillations similar to the B-C oscillations 
shown in Fig. [5] were also reported in this experiment. 
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Therefore, our theoretical results give an alternative ex- 
planation for oscillations in electronic response observed 
for SEs on superfluid helium 4 He. 

CONCLUSIONS 

In summary, we have analyzed the nonlinear WS trans- 
port over superfluid 3 He and 4 He under AC conditions. 
The theory developed for a spatially uniform alternating 
current indicates that the field-velocity relationship ob- 
tained previously in the classical DC model of B-C scat- 
tering is not applicable for time averaged quantities such 
as the first Fourier coefficient. The detailed analysis is 
given for two important limiting cases of low and high fre- 
quency of the electron current. For frequencies which are 
much lower than the ripplon damping coefficient, calcu- 
lations based on the new theory lead to asymmetric B-C 
peaks of a saw-tooth shape which are strongly broadened 
at the right side. The broadening of the right-side tails do 
not depend on small ripplon damping. The left-side tails 
of the B-C resonances become even steeper which pre- 
serves the main B-C anomaly in the field-velocity char- 
acteristic. 

For current frequencies which are comparable with the 
ripplon damping or even higher, the new nonlinear regime 
of B-C scattering of the WS is reported. In this regime 
each B-C peak is transformed into an oscillatory field- 
velocity relationship due to interference of ripplons mul- 
tiply excited at different times. The evolution of surface 
displacements of the dimple sublattice with increasing 
the current amplitude calculated in this work, as well as 
the new field-velocity relationships obtained for alternat- 
ing current, help to understand the nonlinear conductiv- 
ity of the WS on superfluid helium observed in different 
experiments. 
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